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Detailed experimental and theoretical studies of the low-temperature specific heat, magnetic 
susceptibility, thermal expansion and magnetostriction of the orthorhombic compound CeNiSn are 
presented. All anomalies observed in the thermodynamic and magnetic properties of CeNiSn are 
explained in a framework of a model of metallic Kondo lattice with well developed spin-liquid-type 
excitations. The pseudogap behavior of these excitations appears due to interplay between spin- 
fcrmions and soft crystal-field (CF) states. The thermodynamic relations for the spin liquid are 
derived. Together with the explanation of inelastic neutron scattering spectra given earlier within a 
same approach these studies of the low-temperature thermodynamics and magnetic response give a 
consistent description of the nature of anomalies in the low-temperature thermodynamics of perfect 
and imperfect CeNiSn crystals. 

PACS numbers: 75.10.Dg, 75.30.Mb, 71.70.Gm 

I. INTRODUCTION 

The orthorhombic compounds CeNiSn and CeRhSb are known as Kondo lattice systems with peculiar thermo- 
dynamic and magnetic properties. Unusual features are observed at low temperatures T < T* in the specific heat, 
the thermal expansion coefficient, the magnetic susceptibility, the magnetostriction, and the NMR relaxation rate 
(see 0J^] for a review of early data). The characteristic temperature T* is ~ 10K for both systems. It should be 
emphasized that this temperature is much less than the Kondo temperature Tjc estimated by standard methods, e.g., 
extracted from the logarithmic high-temperature dependence of the electrical resistivity. In the early measurements 
the electrical resistivity showed an upturn at low temperatures in the temperature region T < T*. This upturn was 
interpreted as the indication to a non-metallic ground state of these systems, and the energy gap in the heavy electron 
spectrum was claimed to be responsible for the peculiar behavior of CeNiSn and CeRhSb. These materials together 
with the cubic Ce- and U-based compounds, Ce3Sb4Pt3 or UsBi4Pt3 were classified as "Kondo insulators" M. 

Later on it turned out that significant differences exist between the real-gap cubic semiconductors and the or- 
thorhombic CeNiSn family (see [Q for a review). Most striking was the observation that the CeNiSn single crystals 
of good quality show metallic character of the resistivity |5[ at very low temperatures, and such behavior seems to 
be incompatible with the idea of a gap or pseudogap in the electron spectrum. Comparing the metallic behavior of 
electron transport with the anomalous low-temperature thermodynamics, one could suspect that the electronic spec- 
trum with the pseudogap used in early phenomenological theories hardly can be responsible for all low-T peculiarities 
observed in the physical properties of CeNiSn and CeRhSb. Meaningful, e.g., is the fact that the unusual temperature 
dependence of the NMR relaxation rate 1/T\ ~ T 3 which was explained by the V-shape form of the density of electron 
states around the chemical potential at the bottom of the pseudogap Q is observed exactly in the same temperature 
interval where the conventional Fermi-liquid-type T 2 law is seen for the electrical resistivity ||. At T < IK the 
relaxation rate obeys the linear-T Korringa law characteristic for fermions with constant density of states Q. One 
more striking feature of the low-energy excitations in CeNiSn is the extremely complicated (Q, w)-dependent structure 
of the inelastic magnetic scattering spectra that was observed in the same temperature region T < T* Jt],|| . In gross 
features these unusual spectra also can be interpreted in terms of a pseudogap in the spin-excitation spectrum ||, 
although this phenomenology seems to be too simplistic to explain numerous details of the highly anisotropic neutron 
scattering cross section. 

The theoretical approaches to the problem either implement the idea of a Kondo insulator with all its shortcomings, 
or try to offer alternative mechanisms which are based on a metallic type of electron spectra and seek the explanation 
of low-temperature thermodynamics and magnetic response in the unusual properties of the magnetic excitations. In 
the first case the starting point of the theory is the mean-field slave-boson approximation to the Anderson lattice 
The latest version of mean-field hybridization theory M refers to the actual symmetry of f-electron states 
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in the orthorhombic crystal. Since this procedure implies strong coupling of spin and charge degrees of freedom, the 
gap (or pseudogap) in the excitation spectrum necessarily means a semiconductor or semimetallic type of electrical 
resistivity which, apparently, contradicts the experimental data mentioned above. 

An alternative approach was offered in Ref. fll2fl . In this theory new characteristic features with an energy scale 
of T <C Tff appear in the spectrum of the spin excitations due to the interplay between the non-local (spin-liquid) 
excitations characterized by the energy scale of Tr- and the single-site crystal field (CF) excitations with the energy 
Acf < Tk- Within this model the semi-quantitative description of the low-energy specific heat and the thermal 
expansion coefficient was given in Ref. |fl3|j . The CF levels are not seen directly neither in CeNiSn, nor in CeRhSb 
, and this result indicates that these local excitations are " dissolved" in the continuum of low-energy excitations of 
the Kondo lattice. However, the indirect estimate of the magnitude of crystal field created by the Ni ions on the Ce 
site confirmed the validity of the inequality Acf < Tk- 

Basing on the available experimental data related to the structure of magnetic excitations in CeNiSn, the quanti- 
tative theory of interplay between heavy fermions and CF excitations in CeNiSn was offered in Ref. jl5| . The theory 
involves the idea of spin-liquid excitations of resonating valence bond (RVB) type which transform at low temperature 
into well defined fermions with their own dispersion predetermined by dispersion of RKKY exchange in the Brillouin 



zone |16j. As a result of the interplay between these excitations and soft CF states, the spectrum of spin-fermions in 
the low-symmetry lattice of CeNiSn transforms in such a way that the a deep minimum appears in the spin density 
of states (DOS) in the vicinity of the spinon Fermi level. Since the spin-liquid excitations are decoupled from the 
charged Fermi-liquid excitations in the conduction band, the gap in the spin DOS does not imply a corresponding 
gap in the electron DOS, and the system possesses metallic conductivity whereas the spin excitations are responsible 
for the thermal properties. It was shown in Ref. |l5|] that the inelastic transitions between the spinon states in the 
Brillouin zone are responsible for the complicated picture of inelastic magnetic neutron scattering. The successful at- 
tempt of the quantitative description of the magnetic scattering function S^Q, hu), gives us strong arguments in favor 
of the existence of spin-liquid correlations in CeNiSn and related materials. Moreover, the fitting of the theoretical 
spectra to the experimental S(Q, Tiuj) provided us the values of the model parameters. With these data at hand we 



are able to give a q 
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: uantitative description of the temperature dependence of various thermodynamical functions on 

| 13| that the spin-liquid-type excitations give the main contribution to the low-T thermodynamics. 

Thus, a unified description of the low-temperature and low-energy properties of the orthorhombic CeNiSn family 
becomes possible. 

The main purpose of the present paper is to give a detailed experimental and theoretical picture of the low- 
temperature specific heat, magnetic susceptibility, thermal expansion and magnetostriction coefficients in the CeNiSn 
family. This description should be consistent with the picture of magnetic excitations, as given by the inelastic neutron 
scattering experiments. Some of the experimental data for CeNiSn and CeRhSb were published in Refs. (l^Jl^]. The 
first attempts of describing the inelastic magnetic spectra, low-temperature specific heat, and thermal expansion of 
these systems by using the same model p^|,H demonstrated the validity of the spin-liquid description. 



II. EXPERIMENTAL 



The magnetostriction of single-crystalline CeNiSn was measured in magnetic fields up to 8 T at selected temperatures 
of 0.5, 1.4 and 4.3 K. The magnetostriction is defined by A = {L(B) - L(0))/L(0), where L(0) is the length of the 
specimen along a certain crystallographic direction in zero magnetic field. The magnetic field was always applied along 
the orthorhombic a-axis, while A was measured along the field direction (A a ) and perpendicular to the field direction 
along the b- (A&) and the c-axis (A c ). The volume magnetostriction is defined by A.„ = A a + A^ + A c for a fixed field 
direction. The field was applied along the a-axis because this is the easy axis for magnetization: the low-temperature 
susceptibility \a is about a factor two larger than Xb and Xc, and, moreover Xa exhibits a pronounced maximum at 
12 K. 

The experiments were carried out on a Czochralski grown single-crystalline sample. The sample was shaped by 
means of spark erosion into a cube with edges along the principal axes of the orthorhombic unit cell (a x b x c w 
2 x 2 x 2mm 3 ). The magnetostriction was measured using a sensitive parallel-plate capacitance cell machined of 
oxygen free high-conductivity copper. The magnetostriction cell was fixed to the cold plate of a 3 He-insert, which is 
operated with an adsorption pump. The 3 He-insert could be placed in a superconducting solenoid with B max = 8 
T. The magnetostriction was measured by recording the capacitance, while slowly sweeping the field. Temperatures 
were stabilized by regulating on a field- insensitive Ru02 chip-resistor which served as thermometer. 

The experimental results are shown in Fig. la-c, while the coefficients of magnetostriction A^ = L^dL/dB, obtained 
by differentiating the data of Fig. 1 with respect to the field, are shown in Fig. 2a-c. At the highest temperature, 
T = 4.3 K, \ (i — a, b, c) is a monotonous function of the field. The crystal expands in a — b plane and shrinks along 
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c-axis when B || a. The magnetostriction is anomalous in the sense that the curves A^(-B) deviate from the standard 
linear behaviour for paramagnetic systems. At lower temperatures this anomalous behaviour becomes stronger and 
dLi/dB change their signs with increasing field. For instance, at T — 0.5 K, the a — b plane shrinks till ~ 4.5 T 
and the c-axis expands till ~ 7 T. The anomalous behavior is also seen in the volume magnetostriction as a negative 
contribution at low temperatures, although is not very pronounced. 

The magnetostriction data are in good agreement with previous thermal expansion measurements in zero and 
applied magnetic fields (B \\ a) of 4 and 8 T, taken on the same single-crystalline specimen PHI] . Strong anisotropy is 
observed in the linear expansion (pn = L~ 1 (dLi/dT)): the dependence ot c (T) is anomalous with respect to a a ,b{T). 
In magnetic field a sign reversal takes place at low temperatures (T < 3 K at 8 T ): a c (T) becomes negative, 
while a a and ab become positive. The ai(T)-curves show several anomalies, but the coefficient of volume expansion, 
a v — cn a + ai, + a c , is monotonous. Our magnetostriction data are also in excellent agreement with the data reported 
in Ref. in the temperature range 0.1- 4.2 K and field range up to 20 T. 

It is known that the reversible volume magnetostriction is thermodynamically equivalent to the strain dependence 
of the magnetic susceptibility x(B,T) ||(]] 
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(the magnetic susceptibility is defined as x(E,T) = M(B,T)/B, where M(B,T) is the magnetization). Therefore, 
one can extract from the experimental result the field and temperature dependence of the magnetic susceptibility 
volume derivative (dx(B, T)/<91og V)b,t- Introducing the doubly differential magnetostriction coefficient 

K(B,T) =B- 1 (dXv/dB) ViT (2) 

one can express the logarithmic volume derivative of the magnetic susceptibility as 



T\ \'{,(B,T) 
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It is seen from Fig. 3 that the temperature and field dependence of the volume derivative of the magnetic susceptibility 
deviates from the normal temperature and field independent behaviour at low temperatures and in low magnetic fields. 



III. HAMILTONIAN AND THE ENERGY OF THE SPIN LIQUID 



The orthorhombic compounds CeNiSn and CeRhSb are usually classified as Kondo lattices with moderately heavy 
fermion (HF) properties. The basic Hamiltonian which describes the Ce-based HF systems is the Anderson lattice 
Hamiltonian for the Ce 3+ (/ 1 ) ion hybridized with the conduction electrons. In the Kondo lattice limit when the valence 
of the Ce ion is close to integer, one deals with well localized f-electrons for which the inequality V£ fA <C £f — -Er 
is believed to be valid (here is the hybridization matrix element between the f-electron localized on a site i in 
a state |A) = with the energy Er of the f-electron in a crystal field and the partial component of the Bloch 
wave |fcA), v is the row of the irreducible representation T of the crystal point group, ep is the Fermi energy of 
conduction electrons). This hybridization integral is taken in the Coqblin-Cornut (CC) approximation [p2|| which 
represents the Bloch functions by their partial waves c kA , and takes into account only the diagonal in A hybridization 
matrix elements V£ A = (fcA|F|iA). Then the hybridization effects are reduced to exchange-like interaction between 
the localized f-electrons and the conduction electrons with an effective coupling constant 

jf A '(k,k') = Vt A Vi, A ,/(e k -E r ). 

As was shown in |l2| , |l5| , the non-CC hybridization Vl A — (iAH^'|fcA') is of crucial importance for the interplay 
between the one-site crystal-field excitations and the non-local spin liquid excitations (here V is the component of 
the crystal field which has a symmetry lower than that diagonalizing the f-electron energy terms Er- Respectively, 
the non-CC effective exchange constant is introduced as 

Jt A \k,k') = V^*Vl, AI /{e k -E T ). 

In the case of completely suppressed charge fluctuations in the f-channel the sf-exchange can be taken into account 
in the second order approximation, and one comes to the effective RKKY-like Hamiltonian, where the f-electrons 
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are represented only by their spin degrees of freedom described by the spin-fermion operators fi\. When the CF 
excitations are involved, this Hamiltonian acquires the following form (detailed derivation of H s can be found in Refs. 

Ell): 



H' = Hf + H, 



ff{ c ) i rr(™ c ) 

RKKY RKKY 



Here 



ff J = 5> r |iA)(iA| 



(4) 



(5) 



i.A 



describes the Ce(/ 1 )-ions on the lattice sites. 



AA' 



(6) 



corresponds to effective covalent contribution to the one-site CF splitting due to virtual sf-transitions. Here 



AA 



e k -E r ' 



e k -E r 
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The effective exchange interaction mediated by conduction electrons is given by the last two terms in the Hamiltonian 



i#i' 



HrKKY — E y^^»' A /iA/iA'/i'A'/i'A 



(7) 
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and the non-CC interaction is represented by the last term H^^ky which is responsible for the interplay between 
the HF and CF excitations in our model, 



A#A" 



h^ky = E E [^' A " A 'fLhA'fk<f^ + h. 



ii' AA'A" 



(8) 



This is the lowest in V% term among the non-CC indirect exchange interactions which admix the excited CF states 
| A) =| Ei/') to the ground state doublet | A) =| Gv). 

The uniform spin-liquid state in the Heisenberg-like Hamiltonians with antiferromagnetic exchange constant is 
described by the free energy expression 



T = p- X \ £(/3')d0' - /r x <S 
Jo 



(9) 



where (3 1 = fcgT, is the magnetic entropy at T — > oo, and £ is the average value of the Hamiltonian expressed 
via the two-spinon correlator. In the case of the isotropic Heisenberg Hamiltonian this average energy is given by 



£ = E^(|A, 



(10) 



where 



(ii) 



is the non-local operator creating the RVB pair, a stands for the "flavor" (e.g., spin projection in case of pure spin 
states) . After Fourier transformation the average energy of the uniform spin liquid acquires the form 



£ = oEE^w a p a i> 



(12) 
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Here ipu = J2 n exp(— ikR„) is the structure factor for the exchange interaction. 

Usually, in 3D Heisenberg lattices the spin liquid state has higher energy than the AFM state |24j , and the standard 
mean field approach predicts magnetic order at low temperatures. However, the mechanism stabilizing the spin liquid 
state in Kondo lattices was proposed in Refs. [ fL6l|25fl . It was shown within the mean-field approximation that the AFM 
phase can be suppressed by Kondo-type screening, provided X ~ UbTk, and that the spin- liquid state which is not that 
sensitive to Kondo scattering can be realized instead. Recently it was pointed out |g6[ that the influence of low-lying 
relaxation modes in the spin system can transform the phase transition to the spin liquid state into a crossover. The 
low-lying excitation mode (in particular, the soft CF excitations) can play a similar role in stabilization of the spin 
liquid state, and the thermodynamics in this case should be described by the equation generalizing eq. ( |l2] ) for the 
case of CF excitations admixed to the ground Kramers state of the rare-earth ion. The dynamical correlation function 
(A p A q ) w determines the frequency dependence of inelastic magnetic neutron scattering |l5| , so the possibility opens 
for a unified description of low-temperature thermodynamics and the low-energy spin excitations. 

To realize this possibility we first should derive the expression for the free energy of the Kondo lattice described by 
the Hamiltonian H s This means that we should find the energy E or, eventually, to diagonalize the matrix 

M = H' 5 - 1 ■ E (13) 

in terms of the variables A. 

Having in mind the low symmetry of the CeNiSn lattice, we consider the general case of the elementary cell 1 
containing several sublattices £ = l,...,L possessing the same point symmetry group as the Ce ion which total 
magnetic moment is J = 5/2. When diagonalizing (H s ) given by eq. (Q) in terms of spin liquid variables, we 
introduce a single anomalous correlator A G = (/i/j/q/ijg); which corresponds to the ground state doublet A = G. 
In the course of the diagonalization procedure it turns out that this parameter determines the dispersion of both the 
lower and the higher branches of the excitation spectrum which arise due to interplay between CF and HF excitations 
(see Appendix A). 

We introduce the Fourier transformation 

/iU = ^ 1/2 E e4kl0 ^(^ k )/ k , ( 14 ) 

to the basis {k^} which diagonalizes the translationally invariant matrix M (eq. |l^) and the averages (/^/a/i^a) 
(y = 1, (2 J + 1)L). The eigenvectors are orthonormal, 

£e^,k)[ef(^k)]* = <w%, 

V 

£e^,k)[e^,k)r = <w . (15) 

A? 

As a result, the average energy becomes the functional of the "occupation numbers" 

(/L/kv) = «.k!Ak"W , (16) 
and the final equation for energy per Ce ion £(T, {nk}) has the form 

a(0) 

S(T, {n k }) = -jZZ Yl n k ,<MW) (17) 

ki/ 

(see Appendix B). Here Agl is the energy of the lowest CF excitation which is introduced for the sake of convenience 
to make all matrices dimensionless), {nk y } is the set of average occupation numbers which obeys the mean-field 
global constraint condition 

A^^n^l ■ (18) 



ki/ 



The occupation numbers 



n-ki/ 



1 + exp 



knT 



(19) 
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are defined in terms of form factors <I>kiy 

^ = EE e ^' k )#/( k ) K'^' k )]* ■ (20) 

ee aa' 

( fj, is the chemical potential). Then the matrix Z represented by its matrix elements 

Z$ (k) = ^ AA V + i E elkUA #' (»)> ( 21 ) 
should be diagonalized to find the eigenvectors A (£, k). Here u = 1 — 1'. The matrix F AA has the form 

^ AA ' = 6 AA , (ea + BAA + J2 4' 0- j / A cf + S AA '/A^ (22) 

(7^,(1 — 1') is the dimensionless exchange integral, see Appendix A.) Finally, the variables A^A (u) describing the 
RVB state [cf. eq. ([Ti"l)1 are to be obtained self-consistently from the system of equations (|gh|) . 
Thus, to calculate the thermodynamic coefficients of CeNiSn we use the following procedure. 

(i) We find the eigenstates £ki/ of the matrix <I>ki/ which depend on the parameters of the Hamiltonian H s . 

(ii) These eigenstates are used to calculate the imaginary part of the correlation function /C(Q, ui) = (JqJ-q) w 
which determines the dynamic magnetic response of the system, 

Im/C Q/3 (Q, w) = E E "k^ 1 _ n *-Q,») ( k » I Ji I k + Q> A ( k + Q. v ' I J/3 1 kiy ) 5 + ^ - ^k+Q,*') (23) 

and, therefore, scattering function of magnetic neutron scattering [ p^| . 

(iii) The average energy £ is determined by the Fourier component of this correlation function taken at zero moment 

JC { ° ] = J dkdu1C a P(k,Q,uj) (24) 



(sec cq. |l2| and Appendix B). Being diagonalized in terms of the eigenstates |ki/), this energy is given by eq. (|l^ 

(iv) Then, using experimental results of neutron scattering and thermodynamic measurements we fit the model 
parameters to describe the neutron scattering function in absolute units and the heat capacity. Since the heat 
capacity, unlike the neutron scattering spectra, is sample dependent we obtained two set of parameters. The first 
set describes the data for high quality samples and the second one corresponds to the sample used in dilatometric 
measurements. 

The uniform static susceptibility x(^) characterizes the thermodynamic response to an external magnetic field and 
its volume dependence (magnetostriction is determined by the limiting value of /C(0, 0). We calculate it from the 
relation x(T) — M/B where B is magnetic field, and M is the magnetization of the spin liquid. Only the Zeeman 
mechanism of this magnetization is taken into account in the case of weak fields (fJ-sB less than the characteristic 
coupling parameters which determine the spin-fermion spectrum). 



IV. THERMODYNAMIC RELATIONS 



We suppose that the low-temperature thermodynamics of CeNiSn is determined mainly by the spin-liquid excita- 
tions. Since the spin liquid is an unconventional Fermi liquid, and since our treatment of the spin-liquid state inherits 
some of the shortcomings of the mean-field approximation, we start the discussion of the thermodynamic relations 
with a more detailed analysis of spin entropy. It is well known p7[ that one should take special precautions to 
eliminate the unphysical states when introducing the fermionic representation for the spin operators (e.g., the states 
doubly occupied by fermions with opposite spin projections which are absent in original spin representation should be 
excluded). Without such exclusion the wrong temperature behavior of entropy S(T) will result in incorrect descrip- 
tion of the specific heat and other coefficients which are connected with the specific heat by strict thermodynamical 
relations. 

To verify the applicability of our approach we compared the number of states in our model with that in the 
usual Fermi liquid. The model situations considered in Appendix B demonstrate non-universality of the S(T) law 



G 



and its crucial dependence on the parameters of the Hamiltonian, and, in particular, on the character of admixing 
the magnetic CF excitations to the lowest Kramers doublet in the course of forming the spin-fermion branch of 
the excitation spectrum. The diagonalization procedure described above gives the equation for the average energy 
which is sensitive both to the degeneracies of the bare states of the Hamiltonian Hf which are lifted by the spin- 
liquid correlations, and to the temperature compared with the degeneracies lifted already in Hf (CF level splitting). 
According to the calculations for the third model of Appendix B which is close enough to the real situation in CeNiSn 
(see below) the average energy of the spin liquid (|lj]) in the Kondo lattice with several sublattices and several Kramers 
doublets involved can be presented in the form 

A (0) 1 

£{T) = -^f^(T)J2^U{n k }) (25) 



where A.q F = Ee — Eq, and the factor \n(T) reflects the abovementioned fundamental difference between the 
spin-fermion state and the conventional Fermi liquid. 

This factor is essentially nonuniversal: it depends both on the degeneracy of the low-energy branches and the 
lattice geometry. We simplify our consideration by adopting a single value of this parameter for a given geometry of 
the lattice and given set of model Hamiltonian parameters. As is shown in Appendix B, the main quantity which 
predetermines the effective value of k at low temperatures, k^T <C App, is the degeneracy lifted by spin liquid 
correlations. According to the results of description of the neutron scattering spectra in CeNiSn |fl5|| , the "hidden" 
degeneracy of the spectrum equals 4 because (i) all Ce ions are in equivalent crystallographic positions, (ii) the mixing 

interactions comparable in magnitude with £±§p connects only Ce ions belonging to the same fec-planes (see fig. 4) 
although the CeNiSn lattice formally has four sublattices, (iii) only one excited state \E±) = | ± 3/2) interplays with 
the ground state |G±) = a\ ± 1/2) ± b\ T 5/2) which is responsible for the formation of the spin-fermion branch of the 
spectrum. As a result we come to the situation with two Kramers doublets and two sublattices which is treated in 
details in the last example of Appendix B. Therefore we adopt the value of n = 1/4 at k^T <C A[?^. 

At high temperatures the admixing of higher states (the magnetic CF excitons and the branches split due to 
inter-sublattice exchange) becomes essential. As a result, the estimations of the coefficients in equation ( |82j ) together 
with the normalization condition ( p"5| ) give the value of k w 1/2 for these temperatures. Eventually, at high enough 
temperatures exceeding all energy splittings in our Kondo lattice the normal behavior S(T) is restored, but the 
apparent p8[ | entropy deficite at low and intermediate temperatures is an intrinsic property of the model. This 
deficite is an observable effect and its existence was noticed in many measurements (see below). 

The spinon contribution to the heat capacity cy (T) per mole of Ce ion for fixed volume V was calculated using the 
standard expression 

c v (T) = N A (d£(T)/dT) v (26) 

where £ (T) is the spinon energy per magnetic ion ( p5| ) and N A is the Avogadro number. The molar spinon entropy 
for fixed volume S m and the free energy of the mole J- m (T) are found from equations 

S m (T) = [ T °-^-dr , (27) 



T m {T) = N A £{T) - TS m (T) . (28) 

It is known that, in general, such thermodynamic characteristics as the magnetic susceptibility, the volume expansion 
coefficient and the volume magnetostriction can be estimated from the dependence of free energy J-" m (T, V, B) on the 
volume V and magnetic field B. We find these dependences within a framework of the model of spinon spectrum 
which was successfully used in fl8|| for the description of inelastic neutron scattering spectrum. 



V. MODEL OF SPINON SPECTRUM 

CeNiSn crystallizes in the orthorhombic lattice which belongs to the noncentrosymmetric space group Pn2\a [f29| . 
The point symmetry of the crystal field on Ce ions can be treated as nearly trigonal (Dzd) with the rotation axis 
parallel to the a-axis of the crystal, and the monoclinic distortion (C s ) can be considered as a small correction to 
the trigonal crystal field |1J]. Therefore, to describe the bare CF states A we use the irreducible representation of 
the trigonal point group D^d- It is shown by indirect experiments H] and confirmed by the quantitative agreement 
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of the calculated and experimental inelastic neutron scattering spectra pH] that the ground state level and the first 
excited level form a pair of Kramers doublets 

\G±) =o|±l/2>±6|=F5/2) , (29) 
\E±) = |±3/2) (30) 

separated by the energy interval Aqf < 40X which is much less than the energy of the second excited CF level of Ce 
A^p. The recently reported excitations centered around 40 meV |fj(i|] , apparently should be ascribed to this second 
CF state. 

To calculate the average energy £(T) one should solve the system of equations (^,^l|,|2[|9|) under the constraint 
(|l8|). We are interested in the low-T thermodynamics of CeNiSn, (T < 20K). These temperatures are essentially less 
than the bandwidth of the RVB band (W sw 150K, sec |fL5[ ), so we can treat approximately the anomalous averages 
A GG as temperature independent correlators ]3l| , |3^ |. 

It was shown in p5| that the interactions which form the dispersion of the spin liquid excitations are confined 
mainly within the be plane of the CeNiSn lattice. The structure of this plane is determined by the orthorhombic 2D 
elementary cells which contain two Ce ions in the sites i = 1£ where £ = 1, 2 is the sublattice index (see fig. 4). This 
network is defined by the Bravais vectors B = (b, 0) and C = (0, c) and the basis vector d = (0, —6/2, c/2 — O). Here 
O is the orthorhombic distortion which transforms the one-ion hexagonal lattice into the two-ion orthorhombic one. 

To describe the 2D spinon spectrum at low temperatures T — *■ one has to introduce the coupling constants which 
describe the matrix Z^ (k) (|l]) . We confine ourselves with the simplest nearest neighbours (nn) approximation and 
introduce the parameters Tw which describe the in-sublattice 

Ti = A^ F A GG (nn)/2 = A*>« G (m)/2 (31) 

and inter-sublattice 

T 2 = A^A GG (nn)/2 = A^A GG (nn)/2 (32) 

coupling. These constants are responsible for the formation of the spinon spectrum which arise due to RVB correlations 
within the ground state CF level A = G pq] . Similarly, the interplay of spinons with the lowest excited CF state 
A = E is defined by the intra-site nondiagonal matrix element given by eq. ( p2| ) 

G x = F GE = F EG (33) 

and the inter-site mixing coefficient 

G 2 = AP F Ag E (nn)/2 = A^Af G (nn)/2; £ + £' (34) 

The condition £ ^ £' for intermixing of ground and excited doublets arises due to the orthorhombic distortion which 
results in the two-sublattice structure of the plane (Fig. 4). 

To introduce the renormalised CF splitting Acf one has to consider the diagonal terms of eq. (p2[). This contri- 
bution 

A-cf — A^ F (F ee — _F GG ) (35) 

renormalises the bare value of CF splittng A^ = Ee — Eq Then two terms determine the renormalization of CF 
splitting A CF 

Acf = A^p + Sx + S 2 . (36) 

The first correction 

St = B EE - B GG 7 (37) 
is defined by intra-site processes, and the second one is due to inter-site interaction 

= E{ J a'( 1 - 1 ') -J&a -0} ■ ( 38 ) 

Details of the procedure which diagonalizes the matrix (k) and gives the eigenvectors (£, k) and eigenvalues 
$ki/ are described in pa]. 
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VI. SAMPLE DEPENDENCE OF THERMODYNAMIC PROPERTIES 



It is known that the thermodynamic properties of CeNiSn depend on the specimen quality |33|| . At the same time 
the main features of inelastic magnetic scattering spectra are the same (i.e. 2.5 meV and 4 meV inelastic excitations) 
for different specimens [M|[j34[. Therefore, the theory which consider the thermodynamics should explain both the 
sample dependent thermodynamic properties and the sample independent neutron scattering function. 

First, it is known that the more imperfect is the sample of CeNiSn, the higher is the residual resisitivity [0 and the 
lesser is the mean free path of the electrons A. The increase of A in more perfect crystals results in the decrease of 
the exchange RKKY-type interaction ^(^Ce-Ce) ~ ex P( — ^Ce-Ce/^) (HI wn i c h determines the parameters T\ and 
T%. Therefore these parameters should be smaller for less perfect samples. Moreover, since the relative change of 
the resisitivity is higher along the c-axis than along the &-axis, the suppression of the exchange interaction along the 
c-axis should be larger. Hence, the reduction of 72 is greater than that of T\. 

Second, the increase of the defect concentration results in the increase of the inter-site mixing parameter Qi due 
to lowering of the lattice symmetry. The change of the inter-site mixing parameter Q2 is influenced both by increase 
of defect induced mixing and by reduction of inter-site exchange. Therefore, the influence of imperfections on this 



lg,{\-\') » 



the renormalisation of CF 



parameter is not known a priori. Since it is supposed |lq| that 

splitting Acf — ^cf * s defined mainly by the exchange interaction ^/(l — I'). One of the conditions of spin liquid 
RVB state formation is the positive antiferromagnetic sign of these interactions. Therefore, the renormalised value of 
CF Acf is lower for better samples. 

The last effect which has to be considered is the change of the coefficients of the wave functions (^9|) . This change 
is connected with the renormalization of the wave functions of low-symmetry system when the exchange interaction 
changes [ j36f . 

Two sets of parameters which take into account these tendencies for the high quality specimen 0] and the specimen 
which was used in linear expansion and magnetostriction measurements Jf3|,|l7j are presented in Table 1. The fragments 
of DOS for both sets of parameters are shown in Fig. 5. The heat capacity for perfect and imperfect specimens are 
presented in Fig. 6. To illustrate the relative insensibility of neutron spectra to the sample quality we calculated the 
scattering function (in absolute units) for momentum transfers where 2.5 meV (Fig. 7a) and 4 meV (Fig. 7b) are 
observed. It is seen that the neutron scattering spectra for both specimens coincide in gross features. The experimental 
absolute values || of the scattering cross sections are also reproduced. More detailed description of the influence of 
sample quality on the properties of CeNiSn will be presented in a forthcoming publication. 

In the following analysis we use the set of parameters presented in the second column of Table I which correspond 
to the heat capacity of the less perfect specimen used in dilatometric measurements [[13|Jl7|] . 



VII. MAGNETIC SUSCEPTIBILITY 



The standard definition 



B \ dB 



T,N 



was used in the calculation of spinon contribution to molar magnetic susceptibility Xm- Since the magnetic response 
of CeNiSn is maximal for field B a applied along a-direction [|lj , we consider only this easy axis component of Xm and 
limit ourselves by Zeeman mechanism of polaization of spin liquid described by the Hamiltonian 

Hz = 9jy B JaB a . (40) 

Here gj = 6/7 is the Lande factor for Ce 3+ configuration, vb is Bohr magneton and J a is the a-projection of total 
angular momentum operator. It is known pj that the Sommerfeld coefficient 7 sp = C{T)/T is practically constant 
at lowest temperatures T < IK. Therefore we use the Fermi liquid relations for spinon subsystem at T — ► 0. For 
example, the Wilson ratio for spinon liquid can be derived in close analogy with electron Fermi liquid, expression. 
The spinon wave function in the band v with a wave vector k can be represented as a linear combination 



(41) 



where C\(v, k) are coefficients which obey the orthonormality relations. Therefore, applying the Zeeman operator 
( [40| ) to this state we calculate the factor G sp 



9 



G sp =(^2\^M 2 \ 1 j (42) 

which appears in the Wilson ratio instead of the electron ^-factor: 

Xm (T^Q) = \^-G splsp . (43) 

Neglecting in the simplest approximation admixture of the state | E±) =| ±3/2) to the spinons which form the 
RVB band, we find that the lowest spinon state generated by the level ( p9| ) gives 

G SP = gj (5b 2 /2 - a 2 /2) 2 . (44) 

Then, using the parameters from the last column of table 1, we find the value of 1.34-10~ 3 emu/mol for the spinon 
magnetic susceptibility at T — ► 0. This calculated value is significantly lower than the measured one. Therefore, 
one should conclude, that either magnetic susceptibility is determined not only by spinon contribution, or that the 
mechanism of magnetic field influence on the spin liquid can not be reduced to simple Zccman splitting. 

To resolve this alternative, we calculated the temperature dependence of magnetic susceptibility up to T = 20K by 



direct use of eq. (39). It is seen (Fig. 8) that good agreement with experiment can be obtained if one consider the 
total magnetic susceptibility as a sum of spinon Zeeman part X sp {T) and backgroung contribution \b = 6.35 • 10~ 3 
emu/mol which is constant at T < 20K. This calculation perfectly reproduces the position of maximum and the shape 
of the curve at T < 20K. Since % sp (T) in this temperature interval is determined by the low-energy sharp features of 
spinon spectrum one can conclude that the Zeeman splitting of the structured part of spinon spectrum is responsible 
for observed behavior of Jt(T) . The contribution \b is, apparently, detemined by energy scales which are higher than 
the CF splitting. Possible sources of this contributions will be discussed in concluding section 9. 

It should be noted that the spinon response perfectly reproduces the intensities of the 2.5 meV and 4 meV peaks in 
the inelastic magnetic scattering spectra of neutrons in absolute units. Both peaks are connected with the low energy 
structured part of spinon spectrum. This observation gives one more evidence in favor of the assumption that the 
constant contribution to static magnetic susceptibility is connected with larger energy scales. 

More evidences in favor of additional contributions to magnetic response are provided by detailed study of thermal 
expansion and magnetostriction presented in the next section. 

VIII. TEMPERATURE AND FIELD DEPENDENCE OF LATTICE DISTORTION 

A. Thermal expansion 

The conventional phenomenological analysis of the thermal expansion of Kondo lattices is based on the assumption 
that the characteristic temperature Tk scales all thermodynamic quantities at low T, so the main contribution to the 
volume dependence of these quantities may be characterized by a Griineisen parameter jt k = d log Tic/d log V. We 
have seen that the interplay between heavy fermions and crystal field excitations introduces an additional characteristic 
energy scale Acf and a corresponding coupling constant, so this interplay rules out a possibility of being content 
with a single scaling parameter. Moreover, the energy T* is, apparently, not related directly to Tjc in the CeNiSn 
family. Thus we start this section with the derivation of Griineisen parameters which characterize the spin liquid in 
Kondo lattices with soft CF excitations, still confining ourselves to T < Acf <C TA'.The volume thermal expansion 
coefficient 

«v - (W) (45, 



dT 



p 



(P is pressure) can be expressed in terms of the isothermal compressibility kt = (9 log V/8P)t and the isothermal 



derivative of the total entropy Sy with respect to the volume V 37 



Since the Cornut-Coqblin transformation decouples spin and charge degrees of freedom, the total entropy can be 



expressed as a sum Sy 1 — Sy + Sy of the spinon Sy and the conduction electrons Sy contribution. The entropy 
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Sy = V-f el T of conduction electrons is proportional to the Sommerfeld coefficient j el at low T, which results in a 
linear-T law for the thermal expansion 

a v _ . el (i , #log7 



, K ""\ 1+ m^v) < 47 > 

(see, e.g., [[§). 

Due to the inequality Acf <C Tk the spinon component of ay /T can be decomposed into temperature dependent 
and constant terms, ft is convenient to express the spin entropy and its derivative (^6|) in molar units, 



T 



The isothermal compressibility and entropy per mole {KT[mJ / mol] and S m [mJ j '(mol ■ K)\) enter this equation. Since 
the spinon entropy is a function of model constants Vi — Acf, T\^2,Q\-,Qii one can express its volume derivative in 
terms of the corresponding Griineisen parameters. 

tr = P%. (49) 
a log V 

Strictly speaking, the volume dependence of Ac f has to be expressed in terms of magnetoelastic Hamiltonian constants 
p8[ . Hovever, in case of significant contribution of exchange interaction the standard magnetoelastic Hamiltonian has 
to be revised |56j. Therefore, we prefer to describe the magnetoelastic coupling in terms of Griineisen parameters, 
which reflect the main features of magnetoelastic interaction. 

As a result, the spinon contribution to the volume thermal expansion coefficient acquires the form 

a'y _ S m ^ v ( dlogS, 



r-*¥Vn^>ir- (50) 

The peculiar features of a S y /T are determined by those terms which demonstrate appreciable temperature dependence 
of the logarithmic derivatives d log S m /d log Vi ■ The logarithmic derivatives d log S m /d log 71,2) which can be expressed 
in terms of conventional Griineisen parameter ^t k = dlogTx / dlogV , are temperature independent in the considered 
temperature range T <C Tk- Therefore, these terms can be incorporated into a constant Wy together with the 
temperature independent contribution of conduction electrons 

elf, , <91og7 e '\ S m v f d\ogS„ 



dlogV J T \d\ogT K 



^ = ^'1 + ^7 +«r^7k • (51) 



T 



Finally we come to the following expression for the volume thermal expansion coefficient 

a v w . kt / / dS m \ ( dS m \ ( dS m \ \ 

Alhough one can not distinguish between electron and spinon contribution into the constant term Wy, careful 
analysis of temperature dependent contributions provides important infomation about the volume dependence of 
spinons-CF coupling. It is obvious that the on-site mixing strength Q\ is the parameter which is less influenced by 
volume change than the crystal field splitting Aqf and the intersite mixing parameter Qi- Therefore, one can assume 
that 70-l ~ and then analyse the thermal expansion in terms of two Gruneisen parameters, i.e. 7a cf an( i 702 ■ This 
means that even in the temperature range T < Aqf there are two mechanisms of the volume dependence of the 
spinon spectrum. Therefore, one can expect that any attempt to describe the volume-dependent properties of CeNiSn 
by means of a single Gruneisen parameter will result in the temperature dependence of the latter [ fL9| . 

Our two-parameter procedure gives good agreement with experimental data [ fl7| for the set of parameters 7a cf = 
—90 and jg 2 = —60 (see Fig. 9). It should be noted that the signs of both Gruneisen parameters are reasonable, i.e. 
expansion of the lattice leads to softening of Aqf and to decrease of intersite mixing parameter C/2- 

To analyse the anisotropy of thermal expansion one has to introduce the axis-dependent Gruneisen parameters 

x = d log Acf = a & c (53) 

Acf Q log x 3 
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and 



with the obvious property 



<91og£i 

la' = -w\ ; Xj =a,b,c 

yi o log Xj 



(54) 



a,b,c 

Yi = i = ^CF,Gl 



(55) 



A comparison of the calculated linear expansion coefficients with experimental data |l3[ is presented in Fig. 10. 
Since the intcrsitc mixing parameter Q2 can be attributed to the deviation of Ce sublattice symmetry from hexagonal 
one, one can expect that the magnitude of intersite mixing should be proportional to orthorhombic distortion of Ce 
sublattice (see, e.g., Fig. 1 in |lq| ). Therefore, since the lattice expansion along the c-axis is most sensitive to the 
orthorhombic distortion, the inequality | jg 2 \^>\ ~fg I should be valid. Indeed, the best fit is obtained for -fg 2 — — 60 
and 7^ 2 = jg 2 — 0. The same reasoning leads to conclusion about opposite signs of the influence of expansion in the ab 
plane and along the c axis on the CF parameters. Indeed, the values of fitted Gruneisen parameters are 7- = —60, 



7a 



= -84 and 7? = 54. 



B. Magnetostriction 

Two-parameter Gruneisen analysis of thermal expansion can be used as the basis for quantitative explanation of the 
magnetostriction which in our model is the quantity characterising the sensitivity of spinon-CF interplay parameters 
to the volume change. Reversible volume magnetostriction is thermodynamically equivalent to the strain dependence 
of the magnetic susceptibility x(B, T) ||(J 

W)-^X(B,T,(«2) t . m 

Like the thermal expansion, \' v in low magnetic fields can be decomposed in the sum of electronic X' v el and spinon 



contribution Ay 



Ay — Ay ; + Ay . (57) 



The contribution of conducton electrons is temperature independent and linear in magnetic field for fc^T -C £f and 
v B H < e F . 

Aye/ = ktW$B (58) 

where 

(here \p is the Pauli paramagnetic susceptibility). Therefore, to reveal peculiar features of the differential magne- 
tostriction it is convenient to compare experimental and theoretical results for the doubly differential magnetostriction 
coefficient 

Ay = ^ (60) 

which does not depend on temperature and field in the conventional Fermi liquid. This quantity, is proportional 
to the volume derivative of the magnetic susceptibility (provided the field and the temperature dependence of the 
isothermal compressibility is neglected). 

The Gruneisen analysis of the magnetostriction is similar to that of the linear expansion. One should take into 
account only Gruneisen parameters for which the logarithmic derivative of the magnetic susceptibility demonstrates 
sharp temperature and field dependence whereas the structureless contribution due to the standard Gruneisen pa- 
rameter jt k can be taken into account by the renormalization of the normal Fermi-liquid contribution constant Wy 
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— > Wv- However, in addition to the parameters jg 2 and 7a cf one should take into accout the volume dependence of 
the coefficient a in the wave function | G±) ( p9[ ) 

7a = tt, — 77 ■ (61) 
u log V 

Then the final expression for the logarithmic derivative of the magnetic suceptibility acquires the form 

SX -W V + U ,(J^-) + -,J^A ) (62) 



dlogV I ' ACF VaiogAcF/ T \d\o g g 2 J T '"\dloga 

where Wv incorporates all contributions, which are temperature independent at T < 20K 

Wv=W^ lTK {^-) (63) 

The logarithmic derivatives of the magnetic susceptibility are presented in Fig. 11. The fit of the temperature 
dependence of the volume magnetostriction with 7^ CF and j g found above, W — —1590 TO emu/mol, and 7 a = 230 
gives a reasonable agreement with experimental data at low magnetic fields (Fig. 12). It turns out that the last 
term of eq. ( |62"|) dominates in magnetostriction because the contribution of the first two terms gives a value which is 
significantly smaller than the experimental data. 

Like in case of the magnetic susceptibiltythere is a large constant contribution in additon to the temperature 
dependent term, which may be ascribed to the excitations with the energy scale essentially exceeding that of the 
structured part of the spinon spectrum. Moreover, this background, which is present both in magnetic susceptibility 
and magnetostriction, does not appear in linear expansion in zero magnetic field. 

It should be mentioned that the suggested analysis is succesful only at low magnetic fields. Meanwhile, the field 
dependence of the magnetostriction yields evidence in favor of a non-Zecman influence of the magnetic field on the 
spinon spectrum. The experiments in higher fields give almost temperature independent behavior of the logarithmic 
derivative of the magnetic susceptibility (see Fig. 3) whereas the Zeeman term behaves similarly both in low and high 
fields. 



IX. CONCLUSION 



The theory of spin liquid origin of the low temperature anomalies in thermodynamical and magnetic properties 
of CcNiSn and related compounds offered in jl2],|l3| is based on the assumption that the pseudogap features of 
these properties should be ascribed to spin excitations rather than to a non-metallic electron spectrum. Later on, this 
suspicion was confirmed by metallic- like behavior of the resistivity in samples of good enough quality [p| . Nevertheless, 
the spin-liquid theory met the challenge of explaining not only the low-temperature thermodynamics but also the 
fascinatingly complicated picture of inelastic neutron spectra. Explanation of the mechanism of neutron scattering 
offered in [ fl5[ provided us with a set of parameters which determine the spinon spectrum. Thus we came to the 
quantitative picture of a spin liquid which arises as a result of the interplay between intersite resonating valence 
bonds (spinons) and one-site crystal field excitations. In the present paper the experimental data on volume-dependent 
thermodynamical properties of CeNiSn are collected and the quantitative theory of spin liquid is used for interpretation 
of these data. 

To summarize the results of the realization of the above program, one should conclude that the hypothesis that 
the nonlocal spinon pairs determine the free energy of CeNiSn in accordance with eq. (^|) is confirmed by detailed 
quantitative consideration. These pairs, apparently, are responsible both for the structure of low-energy excitations 
which determine the spectrum of inelastic neutron scattering, and for the low-temperature thermodynamics (specific 
heat, thermal expansion). As to magnetic properties (susceptibility, magnetostriction), the real situation is, appar- 
ently, more complicated than the picture given by the mean-field version of spin-liquid theory. We confined ourselves 
by considering the Zeeman polarization of spinon excitations. This mechanism successfully describes the temperature 
dependences x{T) and A'(T) at low enough T where the spinon excitations are still well defined. It fails, however, 
to give the absolute values of these quantities. Besides, the properties of CeNiSn change radically in strong enough 
magnetic field, and this change is also beyond the applicability of our theory. 

At this stage we can only speculate about the presence of another energy scale in the excitation spectrum of the 
spin liquid. To explain the origin of this high-energy scale one should remember that the coherent branch of spinon 
pairs gives only a small contribution to the low-energy part of the spectral density of spin excitations A s (s). In the 
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absence of long range order the dominant contribution to A s (e) is given by the incoherent structurless continuum of 
paramagnetic spin fluctuations. The closest analogy to our case is the structure of spin-fluctuation spectrum in Cu-0 
planes of high-T c materials described by the t — J model at small doping limit |3!^|l|. In that model the coherent 
branch of Zhang- Rice spinon-holon compound particles forms the low-energy part of the spin-polaron continuum. This 
branch gives a small contribution to A s (e) with characteristic energy scale J, whereas the main part of A s {s) is formed 
by the incoherent structureless spin-polaron continuum scaled by the energy < > J. Our crude enough mean-field 
approximation for the spinon spectrum ignores the incoherent continuum. Besides, the influence of magnetic field 
by no means can be reduced to simple polarization of the spin liquid. We hope that the more refined description of 
spin liquid in Kondo lattices will reveal the real role of both external magnetic field and intrinsic antiferromagnetic 
fluctuations in the magnetic response of these systems. 
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XI. APPENDIX A 



To calculate the thermal energy, one should find the thermodynamic average of the sum of Hamiltonians H^, 
H™ KY and H? KKY 



where 



(nc) 



£(T) = £ h (T) + £ C (T) + £ nc {T) 



£h{T) 



A 



(0) 
CF 



NL 



(64) 



(65) 



1{ AA' 



£c(T) 



A 



(0) 



NL 



(66) 



u'se a 



£nc{T) 



A 



(0) / 
CF 



NL 



££(l-5AA0(/l f € AM'A')X 



AA' 



- I'Wlwfw) + [/^ A (l' - 1)] * (fl eA f KA )} . (67) 
Here N is the number of unit cells, primes in the lattice sums mean that the diagonal terms are omitted, 

/£, (1 - 1') = Z$, c , /4 1$ (!-!')= ltSe A ' / A cf , (68) 
and the matrix F AA is defined by the equation ( p2| ) . The quantities A^p (u) are described by the system of equations 

A^'(u) = 2N- 1 Y,n^e-^{5 AA ,I^(u)e^',k) k)] * + 



(i-*aaO (/^ A '(u)e A '(e',k) [e A '(e,k)]* + [/ A ; A (-u)]*e^(e',k) [e A (e,k)]*)} , (69) 

and this system together with equations (|2^), (|2l]) forms the closed set of equations which should be solved self- 
consistently. 
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Although in the general case of low-symmetry lattices with anisotropic exchange constants (u) and JAr (u) , one 
should introduce several variables A^/(u), for the lattices with high enough symmetry one can confine oneself with 
a single parameter A GG which characterises the intersite correlations within the lowest crystal field level. 

For example, in Bravais lattices with P nearest neighbours in equivalent positions the set of parameters should be 
introduced, 

nn 

A AA ' = -Y A AA '(u) . (70) 
z z — ' 

u 

However, if only the lowest crystal field states | G±) are responsible for the "anomalous" intersite correlations described 
by the parameter A GG 

AGG - /G ^E^ k ,|e G (k)| 2 , (7i) 

there is no need in independent nondiagonal variables A GA . All of them can be expressed via A GG by means of the 
factors q A = I AG /I AG < 1, 

A AG = q A A GG . (72) 

Thus, we assume that the symmetry of CeNiSn lattice is high enough to restrict ourself to a single parameter A GG 
which characterises the intersite correlations within the lowest Kramers doublet. 



XII. APPENDIX B 

First, we demonstrate that the diagonalized form ([It]) contains the correct number of states, namely 2N levels for 
the simple case of one-sublattice crystal with spins 1/2 in each site described by the Hamiltonian 

tf S =2££^^'/iV (73) 

ii' vv' 

{y = ± are the spin projections). In this case the average energy of the spin liquid state £ c is given by the equation: 

£c = ^ = Y z £ £ ^ P ~ ^) n P" n ^' ■ ( 74 ) 
ij pq uv< 

Here z is the coordination number for the nn sphere. Using the mean-field definition of the parameter A = (Ay) and 
the spinon energy e p , 

At ^p = £ ^p-q tanh ^ ' e P = /Ae P ' ( 75 ) 
q 

and the property of </? q = 0, equation ( |74|) is reduced to 

f = /^e p n p . (76) 
p 

This means that our problem is thermodynamically equivalent to the problem of spinless fermions, so the limiting 
value of the entropy for this system has the correct value of = AHog2. On the other hand, the naive mean-field 
treatment of the Hamiltonian ( f73| ) results in the effective Hamiltonian 

Hmf =lJ2 £ P» n P» - ^ Nz \ A \ 2 ( 77 ) 
p 

which gives a wrong value of Soo = iVlog4. 
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The nature of this discrepancy is well known. In the spin fermion representation for the spin 1/2, Si = f} u crfi V i 
(where a is the Pauli matrix) the local constraint f} v f% v — 1 forbids simultaneous creation of both up and down 

spin fermions. Since this local constraint is changed for the global constraint N" 1 ^ kj/ / k „./ki/ = 1, one should find 
a procedure which prevents simultaneous creation of " particle" and " hole" in the spinon spectrum when calculating 
the thermodynamic functions to reproduce the correct temperature behavior of entropy. 

More complicated is the situation with the next model example of the Bravais lattice with two CF Kramers doublets 
\Gv), \Ev) and two intersite exchange c oup ling constants I GG and I GE . The system of equations ( |69| ) now describes 
two parameters, A GG (u) given by eq. (16) and 



a ge (u) = lGEj_J2 ^ 1 6 G (k) | 2 = qA 



GG 



(see eq. [72]). Then the matrix Z ( pl| ) has the form 



V 



Gc vk 










\ GG if k 


1 














A GG ^ k 


qA GG 








qA GG <p k 


1 



\ 

) 



(78) 



with normalization condition (15). 

This case can be treated in the same way as the previous one, provided the intermixing of ground and excited states 
is not too strong, i.e., when 



|1 - ±A GG ^ k 



< 1 . 



Then the contribution of the half-filled lowest spin-fermion band to the energy £ is 



i (0) 

*CF 



k v=± 



where 



e k , = ^-A GG tp k ( 1 - 



|l-±A G <V k | 



(79) 



(80) 



(81) 



In this case we also have the compensation of Kramers degeneracy. It should be emphasized, however, that the second 
branch of the excitations generated by the matrix Z (fT8|) is, in fact, the usual magnetic CF exciton band modified 
by the interaction with the spin liquid branch, and its contribution to the entropy can be treated in conventional 
manner, at least at fegT <C ^cp- 

These two examples demonstrate that there is no universal recipe for calculating the entropy in the systems with 
strong interplay between the non-local spin-liquid excitations and the one-site CF excitations. The third instructive 
example demonstrates the importance of accurate treatment of all degeneracies which could be lifted by the spin-liquid 
correlations. Here we consider the two-sublattice crystal with the crystal field resulting in two equivalent Kramers 
doublets for Ce ion in each sublattice. Then the exchange interactions in H c and H nc terms of the Hamiltonian (Q) 
are described by four parameters and if® , where ££' = 1, 2. Then the matrix Z acquires the form 



A GG, 



( A GG </? k i\Y2<Pk 
A GG cp' k A GG <p k 

qA GG ip k q'A GG ip' k 

\q'A GG ^ qA GG <p k 



GG 

q'A GG V ' k qA GG Lp k 
2 
2 



GASfVk g'Af 2 G ^ k \ 



/ 



(this matrix represents one of two Kramers subspaces in the block-diagonal matrix Z = Z + ® Z~). We assume that 
the in-sublattice structure factor cp k and the in-sublattice coupling constants I\\ are the same for both sublattices. 
Here relation between non-diagonal and diagonal elements of in-sublattice and inter-sublattice coupling constants are 
given by q — I G ^ // G and q' = I G 2 E /I G 2 , respectively. In general case the inter-sublattice structure factor ^ yj k . 
Under these assumptions An = A22 and A12 = A21, and two independent spinon parameters are determined by the 
equations 
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A?i G = ^$> k n k ,^ |H!;il.k)| 



A ^ G = 4z E vUk.Jf 2 e?(2, k) [e?(i, k)] * . (82) 



Nz 

In the case of strong intermixing (I^*/ f» 1 and q, q 1 « 1) all the coefficients |Oj?(£, k)| are of the same order, and each 
of them can be estimated as 1/2 y2. As a result, the contribution of the lowest branches at low temperatures 
/csT <C Ajjp can be approximately represented as 

^^A^A^EE'^L (83) 

k v=± 

because in this case all terms in the matrix Z related to lowest branches <; contain the factors A , as was explicitely 



demonstrated in the above example (see eq. 81) 



Again we see that the diagonalization procedure (17) results in a non- universal form of the average energy for 
spin liquid excitations in comparison with corresponding equations for the conventional Fermi liquids, at least at low 
temperatures. 
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Parameters 


High quality 


B.lNOlten et.al. 


Ti 


16.15 K 


14.40 K 


T 2 


24.65 K 


12.40 K 


Si 


1.83 K 


3.00 K 


02 


-6.60 K 


-5.41 K 


A C F 


-2.0 K 


8.0 K 


a 


0.53 


0.62 



Tabic 1. Parameters which define the spinon spectrum in the samples of higher and lower quality 
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Figure captions 



Fig.l Magnetostriction A of single-crystalline CeNiSn for a field directed along the a-axis (bold lines) and elongation 
(or contraction) along the a, b and c-axis (thin solid, dotted and dashed lines, respectively) at temperatures of 0.5 K 
(a), 1.4 K (b), and 4.3 K (c). 

Fig. 2 Coefficient of magnetostriction A' of single-crystalline CeNiSn for a field directed along the a-axis. All notations 
are the same as in fig. 1 

Fig. 3: Field and temperature dependence of the logariphmic volume derivative of magnetic susceptibility evaluated 
from eq. (j|) (the isothermal compressibility is kt — 1.8 * 10~ n m 2 /N according to pl| ): (a) field dependence for 
T = 0.5 K (squares), T = 1.4 K (diamonds) and T — 4.3 K (triangles); (b) temperature dependence for different 
magnetic fields. 

Fig. 4: Orthorhombic CeNiSn lattice, structure of frc-plane. Two Ce sublattices are denoted by black and grey circles, 
respectively. The orthorhombic distortion O (solid arrows) transforms the simple hexagonal lattice into two-sublattice 
orthorhombic one. The in-sublattice interaction 7[ is denoted by solid double arrow. The inter-sublattice interaction 
72 and Q2 is denoted by dashed double arrow. 

Fig. 5: Fragments of the density of states (DOS) of spinon spectra for high quality sample (bold line) and imperfect 
sample (thin line). 

Fig. 6: Temperature dependence of heat capacity. Dashed line represents the spin-fermion contribution j sp , for im- 
perfect sample, and solid line gives the Sommerfeld coefficient with additional contribution from conduction electrons 
Icond = 8 mJ/mol K 2 . The spinon part calculated for high quality sample is presented in the insert. Experimental 
points for perfect sample are taken from || (triangles) and those for imperfect sample are taken from |17| (squares). 

Fig. 7: Scattering functions of inelastic magnetic neutron scattering in absolute units calculated for high quality (bold 
line) and less perfect sample (thin line). 

Fig. 8: Calculated magnetic susceptibility (line) compared with experimental data for specimen # 4 from Ref. |J. 

Fig. 9: Calculated linear coefficient of volume expansion (line) compared with experimental data JlTj (triangles). In- 
sert: temperature dependence of logarithmic derivatives of entropy. 

Fig. 10: Calculated coefficients of linear expansion (lines) compared with experimental data p"3| ] (points) along a- 
(solid line, squares), b- (dashed line, triangles) and c-axis (dotted line, dimonds). 

Fig.l 1: Components of volume derivative of magnetic susceptibility calculated from eq. ( |(32| ) 

Fig. 12: Logarithmic derivative of magnetic susceptibilty calculated from eq. ( |(32| ) (solid line) and extracted from 
experiment |21| (circles) . 
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